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Abstract We discuss the properties of cosmic ray nuclearites, from the point 
of view of their search with large nuclear track detector arrays exposed at dif- 
ferent altitudes, in particular with the SLIM experiment at the Chacaltaya high 
■ altitude lab (5290 m a.s.L). We present calculations concerning their propaga- 

tion in the Earth atmosphere and discuss their possible detection with CR39 
and Makrofol nuclear track detectors. 
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Strange Quark Matter (SQM) could be the ground state of quantum chromod- 
inamics 1 . The initial hypothesis assumed that SQM is made of u, d and s 
quarks in nearly equal proportions (with some electron component in weak equi- 
librium). Recently it was shown that the so called "color-flavor locked" (CFL) 
SQM, characterized by a Cooper-like pairing between different quarks, could be 
even more stable 

SQM is expected to have a density slightly larger than ordinary nuclear 
matter ^H]; the relation between the mass M of SQM lumps and their baryonic 
number A ( ~ one third of the number of constituent quarks) would be 

M(GeV) < 0.93A (1) 

It was hypothesized that "nuggets" of SQM, with masses from those of heavy 
nuclei to macroscopic values, produced in the Early Universe or in violent as- 
trophysical processes, could be present in the cosmic radiation (the so-called 
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nuclearites) |3| 1 . SQM should have a relatively small positive electric charge, 
eventually neutralized by an electron cloud. If the size of the SQM is large 
(corresponding to masses M > 10 7 GeV), some of the electrons could be in 
chemical equilibrium inside the quark core 0]. Nuclearites larger than 1 A, 
(M > 8.4 x 10 14 GeV) would contain all electrons inside the quark core and thus 
would be completely neutral [3J. A qualitative picture of nuclearites may be 
found in 0. 

An upper limit for the flux of nuclearites may be obtained assuming that 
they represent the main contribution to the local Dark Matter (DM) density, 
Pdm 10- 24 g cm" 3 0, 

, _ Pdmv r „s 
max ~ 2ttM ' 1 } 

where v and M are the nuclearite average velocity and mass, respectively. 

The aim of this note is to discuss the possibility to detect nuclearites using 
large area Nuclear Track Detectors (NTDs) at mountain altitude (in particular 
with the SLIM experiment |S1[7||H])- We classify nuclearites in three different 
mass ranges (and sizes). 

• Low mass nuclearites (LMNs), or strangelets, with masses between 
those of ordinary nuclei (A < 300) and a multi-TeV mass. The upper 
bound is an ad-hoc one; an indirect definition of this category could be 
that many of the properties of LMNs would contradict the assumptions 
made in |3] summarized in Section 3. We shall discuss the case of LMNs 
in Section 2. 

• Intermediate mass nuclearites (IMNs) or, simply nuclearites, 

with masses large enough to be well described by the hypothesis made 
in 0, but smaller than about 10 22 GeV (for M > 10 22 GeV nuclearites 
would traverse the entire Earth) . The mass lower limit of IMNs may be 
about 10 s GeV, above which nuclearites could be detected by experiments 
performed in the upper atmosphere (see also Fig. [^J. Assuming that they 
would travel in space with typical galactic velocities ((3 = v/c ~ 10~ 3 ), 
they would be stopped by the Earth, so they would reach detectors only 
from above. The main properties of IMNs are summarized in Section 3. 

• "Macroscopic" nuclearites, with masses M > 10 22 GeV; assuming 
galactic velocities, such nuclearites would traverse the Earth. They differ 
from IMN's only by size. We shall not discuss this case, since the expected 
sensitivity of SLIM (and of similar experiments) will not compete with the 
limit obtained by MACRO in this mass range P ITUIITT] . 



lr The attention of the authors was mostly focused on relatively large mass nuggets of 
SQM, so that the microscopic properties of the nuclearites would not be relevant. This is 
suggested also by the name they proposed for the newly postulated objects, "nuclearites": 
a combination between "nuclei" and "meteorites" . Note that a nuclearite is an electrically 
neutral state composed of a SQM "nucleus" and electrons. 
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Much heavier nuclearites (M > 10 28 GeV) could be observed as abnormal 
sismic events |2] 2 • 

Calculations describing the production (through binary strange stars tidal 
disruption) and the galactic propagation of cosmic ray nuclearites were recently 
published ^7]. The results could be valid as orders of magnitude for the entire 
mass range of interest; we shall use the predicted fluxes (at the Earth level) as 
reference values. 

Searches for nuclearites (mostly IMNs) were performed by different experi- 
ments ^1^1- The best flux upper limit was set by the MACRO experiment: for 
nuclearites with (3 ~ 10~ 3 , the 90% C.L. upper limit is at the level of 2 x 10~ 16 
cm^sr^s" 1 in the mass range 10 14 GeV < M < 10 22 GeV P HUIITT] 3 . 

SLIM is a large area experiment (440 m 2 ) installed at the Chacaltaya high 
altitude laboratory since 2001; an additional 100 m 2 were installed at Koksil, 
Pakistan, since 2003 4 . With an average exposure time of about 4 years, SLIM 
would be sensitive to a flux of downgoing exotic particles (magnetic monopoles, 
nuclearites and Q-balls) at a level of 10 -15 cm _2 sr _1 s _:l . 

2 Low mass nuclearites (LMNs) 

SQM should be stable for all masses larger than about 300 GeV [3]. Nuclearites 
with masses up to the TeV region could be ionized and could be accelerated to 
relativistic velocities by the same astrophysical mechanisms of normal nuclei of 
the primary cosmic radiation (CR). 

LMNs would interact with detectors (in particular NTDs) in ways similar to 
heavy ions, but with different Z/A. There are different calculations in relative 
agreement with a possible candidate with M ~ 370 GeV and a charge Z ~ 14 

M 

In ref. 0] SQM is described in analogy with the liquid-drop model of normal 
nuclei; the obtained charge versus mass relation is shown in Fig. [2 by the solid 
line, labeled "(1)". Other authors found different relations: Z ~ 0.1^4 for 
A < 700 and Z ~ 8 A 1 ^ for larger masses |T5]: this charge to mass relation 
is shown in Fig. ^ as the dashed line, labeled "(2)" . In [2] it was assumed 
that quarks with different color and flavor quantum numbers form Cooper pairs 
inside the SQM (the so-called color-flavor locked phase), increasing the stability 
of the strangelets. In this case, the charge relation would be Z ~ 0.3A 2 / 3 , shown 
as the dash-dotted line in Fig. ^labeled "(3)". 

Several CR experiments reported candidate events that would suggest anoma- 
lously low charge to mass (Z/A) ratios, which could correspond to those ex- 
pected for SQM 4 . Such candidates are reviewed in [2U||^. As strangelets 
with masses not much higher than those of ordinary nuclei could have the same 

2 In ref. 1131 there was a claim of observing a candidate, but it was discarded in 1141 because 
of timing uncertainities. 

3 This is twice the flux limit obtained for relativistic GUT magnetic monopoles 1121 . as it 
refers only to down-going nuclearites. 

4 The calculations presented in this report refer to the Chacaltaya location only. 
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Figure 1 : Low mass nuclearites (LMN) charge versus mass for different hypothe- 
ses discussed in refs. 0JE3I2]- See text for details. 



origin as CR heavy nuclei, their abundances in the cosmic radiation could fol- 
low the same mass dependence, $ oc Af~ 7 5 , [201 - The existing candidates do 
not contradict such an hypothesis. The solid line in Fig. [21 is the expected 
flux versus nuclearite mass, assuming that the above assumptions are correct. 
Obviously, as the mass becomes larger, the production mechanisms for normal 
cosmic rays cannot anymore apply to nuclearites. 

Different nuclearite flux estimates were recently published ^Zj. They are 
based on the hypothesis that large nuclearites (with masses 10~ 5 — 10~ 2 so- 
lar masses) are produced in binary strange stars systems, before their gravita- 
tional collapse. The propagation inside the galaxy considers also the escape, 
spallation (through which smaller nuclearites are produced) and re-acceleration 
mechanisms. Nuclearite decays are not considered, as SQM is supposed to be 
absolutely stable. The predicted strangelet fluxes around the Earth are pre- 
sented in Fig. [21 for "normal" and CFL strangelets as the dashed and the dot- 
dashed lines, respectively. The small differences (that vanish for larger masses, 
when nuclearites become completely neutral) originate in the slightly different 
charge-to- mass ratios (see Fig. ^1 . 

LMNs would be similar to normal nuclei, except their low Z/A values and, 
most likely, different abundances. If they interact with the Earth's atmosphere 
in the same way as CR nuclei, they would not reach experiments at mountain 
altitude. 

Two different (and opposite) theoretical scenarios, both consequences of the 
hypothesis that SQM is more stable than ordinary nuclear matter, were intro- 
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Figure 2: Expected fluxes for LM strangelets in the CR near the Earth. The 
solid line corresponds to the assumption that their abundances follow the same 
rule as heavy CR nuclei (20] • The dashed and dot-dashed lines are from Ref. 
|17| . and refer to "normal" and CFL strangelets. 



duced in order to allow deep penetration of small nuclearites in the atmosphere; 
none of those mechanisms would allow them anyway to reach sea level; such 
objects could be found only in high altitude experiments. 

2.1 Mass and size decrease of nuclearites during propaga- 
tion 

In |20| it was assumed that small nuclearites could penetrate the atmosphere if 
their size and mass are reduced through successive interactions with the atomic 
air nuclei. The proposed scenario is based on the spectator-participant picture. 
Two interaction models are considered: quark-quark (called "standard"), and 
collective (called "tube- like" ) . At each interaction the nuclearite mass is reduced 
by about the mass of a Nitrogen nucleus (in the "standard" model) , or by more 
(in the "tube- like model"), while the spectator quarks form a lighter nuclearite 
that continues its flight with essentially the same velocity as the initial one. Once 
a critical mass is reached (A ~ 300 - 400) neutrons would start to evaporate 
from strangelets; for A < 230 the SQM would become unstable and decay into 
normal matter. In ref. |22I an estimate was made of the sensitivity of the SLIM 
experiment [B]: the mass number of a nuclearite penetrating the atmosphere 
down to the Chacaltaya lab would be one seventh of that it had at the top of 
the atmosphere. 

The CR39 used in SLIM is sensitive to particles with a Restricted Energy 
Loss (REL) larger than 200 MeV g _1 cm 2 5 . Table 1 presents the minimum 

5 The threshold of a NTD depends on the etching conditions. A relatively high threshold for 
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detectable mass for strangelets reaching SLIM. 



A/Z 
hypothesis 


A 

(at SLIM) 


^0 

(top of 
the atmosphere) 


$ 

cm _2 s _1 sr _1 
(Ref. ED 


$ 

cm _2 s _1 sr _1 
(heavy nuclei) 


Liquid-drop 
or CFL 

Rets, mm 


~ 587 


~ 4109 


6 x 10~ 12 {CFL) 
3 x 10" 12 (normal) 


2.4 x 10" 15 


"Normal" 
nuclearites 
Refs. [HIP] 


~ 210 


~ 1470 


9.5 x 10" 12 


2.8 x 10~ 12 



Table 1: The minimum strangelet masses detectable in SLIM assuming 
different A/Z relations. The masses at the top of the atmosphere are 
estimated as in Ref. j^U) > and the expected fluxes are computed as in Ref. |17|. 
and assuming the same mass abundances of ordinary cosmic ray nuclei. 



2.2 Accretion of neutrons and protons during propagation 

A completely different propagation scenario was proposed in |21| . The authors assume 
that small mass nuclearites would pick-up nuclear matter during interactions with air 
nuclei, rather than loosing mass. After each interaction, the nuclearite mass would 
increase by about the atomic mass of Nitrogen, with a corresponding slight reduction of 
velocity. As the mass grows larger, the loss in velocity becomes smaller. They estimate 
that a strangelet of an initial A ~ 64 and an electric charge of about +2 could arrive 
at about 3600 m a.s.l. with A ~ 340 (3600 m is the altitude of a proposed NTD 
experiment in Sandakphu, India |21|L This mechanism would also imply an increase 
of the electric charge of the strangelet, thus an increase of the Coulomb barrier; this 
may be the main difficulty of this scenario. The flux according to |17| would be of the 
order of 10~ 9 cm _2 s _1 sr _1 , and higher. 

About 171 m 2 of the SLIM modules exposed for an average time of 3.5 years were 
removed, processed and analized. No cadidate survived. The 90% C.L. flux upper limit 
for downgoing nuclearites (LMNs and IMNs) is at the level of 4 x 10 -15 cm _2 sr s . 
This would disfavor the "accretion scenario", and most of the hypotheses quoted in 
Table 1. 



3 Intermediate mass nuclearites (IMNs) 

In |3] was postulated that elastic collisions with the atoms and molecules of the tra- 
versed medium are the only relevant energy loss mechanism of non-relativistic nucle- 
arites with large masses, 

dE 2 

^ = ~ apV ' (3) 

CR39 was chosen in order to reduce the background due to recoil tracks, neutron interactions 
and the ambient radon radioactivity. 
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where p is the density of the traversed medium, v is the nuclearite velocity and a is 
its cross section: 



_ ( 7r(3M/47rpjv) 2/3 for M > 8.4 x 10 14 GeV (corresponding to R N ~ 1 A) 
\ it x 10~ 16 cm 2 for lower mass nuclearites 

(4) 

with pjv = 3.6 x 10 16 g cm 3 . As the chemical potential of the s quarks in SQM is 
slightly larger than for u and d quarks, SQM is always positively charged thus the 
cross section for nuclearites with M < 8.4 x 10 14 GeV is determined by their electronic 
cloud. 

The following calculations apply to nuclearites of mass M much larger than typical 
nuclear masses and /? ~ 10 -3 ; effects due to possible ionization or mass variations 
during their flight in the atmosphere are negligible. We also neglect the gravitational 
acceleration of nuclearites by the Earth. 6 

A nuclearite of mass M entering the atmosphere with an initial velocity vo << c, 
after crossing a depth L will be slowed down to 

v(L)=v e-^ J ° pdx (5) 

where p is the air density at different depths, and a is the interaction cross section of 
Eq. [1 

We consider the parametrization of the standard atmosphere from I23| : 

_h H-L 

p(h) = ae b — ae b , (6) 

where the constants are a = 1.2 x 10~ 3 g cm~ 3 and b ~ 8.57 x 10 s cm; H is the total 
height of the atmosphere ( ~ 50 km). The integral in Eq. 3 may be solved analytically: 

_H_ I H-h \ 

pdx = abe <> I e <> - 1 J . (7) 

Fig. |3 shows the velocity with which nuclearites of different masses reach heights 
corresponding to typical balloon experiments (for instance CAKE, 40 km \'24\). possible 
experiments using civilian airplanes (11 km), the Chacaltaya lab (SLIM, 5.29 km 
and at sea level. A computation valid for MACRO |10| (at a depth of 3400 mwe) is 
also included. The velocity thresholds for detection in CR39 (corresponding to REL 
= 200 MeV g _1 cm 2 ) 7 and in Makrofol (REL = 2500 MeV g" 1 cm 2 ) are shown as the 
dashed curves. 

The decrease of the velocity thresholds for nuclearite masses larger than 8.4 x 10 14 
GeV is due to the change in the nuclearite cross section, according Eq. 2] 

An experiment at the Chacaltaya altitude lowers the minimum detectable nucle- 
arite mass by a factor of about 2 with respect to an experiment performed at sea level. 
If the mass abundance of nuclearites decreases strongly with increasing mass this could 
yield an important increase in sensitivity. 

6 Assuming a nuclearite mass of 1 ng (about 5.6xl0 14 GeV) arriving at an altitude of 5000 
m with j3 = 10 — 3 , the gravitational energy gain would represent less than about 1.5 X 10 — 3 
of the energy loss in the atmosphere; for /3 = 10 -4 (near the Makrofol threshold) this ratio is 
about 0.15. 

7 We recall that in the low background conditions of the Gran Sasso Lab, in the case of the 
MACRO experiment the CR39 detection threshold was set to 50 MeV g —1 cm 2 
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Figure 3: Solid lines: arrival velocities of IMNs at different depths versus nu- 
clearite mass, assuming an initial velocity outside the atmosphere of (3 = 1CP 3 . 
The nuclearites are supposed to come from above, close to the vertical direction. 
The dashed lines show the detection thresholds in CR39 (in the SLIM etching 
conditions) and Makrofol. 



The nuclearite detection conditions in CR39 (expressed as the minimum entry 
velocity at the top of the atmosphere versus the nuclearite mass) for different exper- 
imental locations is shown in Fig. |1] In this case, the constraint is that nuclearites 
have the minimum velocity at the detector level in order to produce a track; we remind 
that for all experiments the REL threshold for detection in CR39 is set to 200 MeV 
g" 1 cm 2 (it was 50 MeV g" 1 cm 2 in the case of MACRO). 

Fig. Q>]shows the same for the Makrofol track etch detector. The detection condi- 
tion corresponding to CR39 at balloon altitude (CAKE) is also shown. 



4 Conclusions 

SLIM is a large area NTD experiment, taking data at the Chacaltaya Cosmic Ray 
Laboratory. In this note we investigated the possibility to search for nuclearites with 
SLIM, assuming two nuclearite mass regions. 

Low mass nuclearites could reach mountain altitudes assuming some peculiar 
interaction mechanisms in the atmosphere. They would produce in NTDs tracks sim- 
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Figure 4: Nuclearite detection conditions in CR39, for experiments located at 
different altitudes. 



ilar to those expected from fast monopoles or relativistic heavy nuclei 8 JBj- SLIM 
will reach a sensitivity at the level of about 10~ 15 cm _2 s _1 sr _1 for a flux of nucle- 
arites coming from above. In the absence of a LMN candidate, SLIM would rule out 
the propagation mechanisms hypothesized. SLIM will be also sensitive to different 
strangelet structure hypotheses ("normal" or CFL, different Z/A predictions), and 
could validate or not the production and propagation model proposed in |17|. 

Intermediate mass nuclearites, entering the Earth atmosphere with typical 
galactic velocities might be detected by large area NTDs as SLIM. The minimum de- 
tectable nuclearite mass is very sensitive to the experiment location: an underground 
experiment like MACRO could search for nuclearites with M > 10 14 GeV, the mini- 
mum mass accessible to ANTARES I25| is of few 10 13 GeV. Detectors at ground level, 
or better at mountain altitudes like SLIM, would decrease the mass threshold to few 
10 10 GeV; balloon experiments are needed to reach few 10 s GeV, while lower mass 
searches have to be done outside the Earth atmosphere. 

SLIM is sensitive to non relativistic (/3 < 10 -3 ) IMNs with masses larger than 
3 x 10 10 GeV; the large REL of IMNs in NTDs and their property to produce identical 
tracks in all the detector sheets in a stack could yield experimental signatures with 
low background. 

Nuclearites with masses between few TeV and about 10 8 GeV are not considered 

8 Note that relativistic nuclei present in the CR cannot penetrate the Earth's atmosfere till 
the Chacaltaya level. 
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Figure 5: Nuclearite detection conditions in Makrofol, for experiments located 
at different altitudes. The dashed curve (shown for comparison) corresponds to 
CR39 in CAKE. 



at this time. In the low part of this mass range, they could still be accelerated to 
relativistic velocities by cosmic fields, and could be detected as LMNs in high altitude 
experiments. For larger nuclearite masses, one would expect the velocity to get close 
to the galactic one, and thus such nuclearites would not be able to reach the detec- 
tor. A search for nuclearites in this mass range could, in principle, be done in space 
experiments, such as the AMS detector on board of the International Space Station 

Fig. |S] summarizes the conclusion, showing the accessible nuclearite mass regions 
for different experiments. 
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